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ABSTRACT 

Extended UltraViolet (xuv) disks have been found in a substantial fraction of 
late- type -SO, spiral and irregular- galaxies. Similarly, most late- type spirals have an 
extended gas disk, observable in the 21cm radio line (Hi). The morphology of galaxies 
can be quantified well using a series of scale-invariant parameters; Concentration- 
Asymmetry-Smoothness (CAS), Gini, M20, and Gm parameters. In this series of pa- 
pers, we apply these to Hi column density maps to identify mergers and interactions, 
lopsidedness and now xuv disks. 

In this paper, we compare the quantified morphology and effective radius (i^so) 
of the Westerbork observations of neutral Hydrogen in Irregular and Spiral galaxies 
Project (WHISP) Hi maps to those of far-and near-ultraviolet images obtained with 
GALEX, to explore how close the morphology and scales of H I and UV in these disks 
correlate. We find that xuv disks do not stand out by their effective radii in UV or 
Hi. However, the concentration index in FUV appears to select some xuv disks. And 
known XUV disks can be identified via a criterion using Asymmetry and M20; 80% of 
xuv disks are included but with 55% contamination. This translates into 61 candidate 
xuv disk out of our 266 galaxies, -23%- consistent with previous findings. Otherwise, 
the UV and H I morphology parameters do not appear closely related. 

Our motivation is to identify xuv disks and their origin. We consider three sce- 
narios; tidal features from major mergers, the typical extended Hi disk is a photo- 
dissociation product of the xuv regions and both Hi and UV features originate in 
cold flows fueling the main galaxy. 

We deflne extended H I and UV disks based on their concentration {Chi > 5 and 
Cfuv > 4 respectively), but that these two subsamples never overlap in the WHISP 
sample. This appears to discount a simple photo-dissociation origin of the outer Hi 
disk. 

Previously, we identifled the morphology space occupied by ongoing major merg- 
ers. Known XUV disks rarely reside in the merger dominated part of H I morphology 
space but those that do are Type 1. Exceptions, xuv disks in ongoing mergers, are the 
previously identified UGC 4862 and UGC 7081, 7651, and 7853. This suggests cold 
fiows as the origin for the xuv complexes and their surrounding Hi structures. 

Key words: 



1 INTRODUCTION 

Interest in the outskirts of spiral galaxy disks has increased 
over recent years as these regions are the site of the most 
recent acquisition of gas for these systems (e.g. Sancisi et al. 



* E-mail: benne.holwerda@esa.int 



2008), as well as low- level star-formation (e.g., Dong et al. 
2008; Bigiel et al. 2010a; Alberts et al. 2011, for recent re- 
sults), making these faint outskirts the interface between the 
island universes -the galaxies themselves- and the cosmic 
web of primordial gas. 

The low-level star-formation was first discovered in 
Ha emission by Ferguson et al. (1998) and Lelievre & 
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Roy (2000). After the launch of the Galaxy Evolution EX- 
plorer (galex, Martin et al. 2005), initial anecdotal evidence 
pointed to ultraviolet disks of spiral galaxies extending much 
beyond their optical radius (Thilker et al. 2005b, a; Gil de 
Paz et al. 2005, 2007b; Zaritsky & Christlein 2007). Sub- 
sequent structural searches for these Extended Ultraviolet 
(xuv) Disks by Thilker et al. (2007) and Lemonias et al. 
(2011) find that some 20-30% of spirals posses an xuv disk 
and 40% of SOs (Moffett et al. 201 1) , making this type of disk 
common but not typical for spiral and SO galaxies. These 
XUV disk complexes are generally ^ 100 Myr old, explaining 
why most lack Ha (Alberts et al. 2011), as opposed to a 
top-light IMF (as proposed by Meurer et al. 2009), and sub- 
solar but not excessively low metallicities (Gil de Paz et al. 
2007b; Bresolin et al. 2009; Werk et al. 2010, 0.1-1 Zq, 
based on emission lines). Additionally, it has been known 
for some time now that atomic hydrogen (Hi) as observed 
by the 21cm fine structure line also extends well beyond the 
optical disk of spiral galaxies (e.g., Begeman 1989; Meurer 
et al. 1996, 1998; Swaters et al. 2002a; Noordermeer et al. 
2005a; Walter et al. 2008; Boomsma et al. 2008; Elson et al. 
2011; Heald et al. 201 lb, a; Zschaechner et al. 2011). In those 
few cases where both high-quality H I and deep GALEX data 
are available, a close relation in their respective morphol- 
ogy was remarked upon (Bigiel et al. 2010a). While we have 
to wait for the all-sky surveys in Hi to catch up to the 
coverage of the galex surveys (e.g., the wide survey with 
WSRT/APERTIF or the WALLABY survey with ASKAP), 
we can compare the morphology for those galaxies for which 
uniform Hi information is available. 

The canonical view of the origin of the xuv disks, is 
that the Kennicutt-Schmidt law (Kennicutt 1998) needs to 
be extended to low global surface densities of gas and the 
formation of individual 0-stars in the very outskirts of disks 
(Cuillandre et al. 2001; Bigiel et al. 2010b). The recent ac- 
cretion of cold gas flows (Keres et al. 2005) into the H i disk 
is the origin for the young stars (the fueling rate implied 
by XUV disks is explored in Lemonias et al. 2011) and Hi 
warps(Roskar et al. 2010). The fraction of spirals that have 
an XUV disk (~ 20 — 30%) supports this scenario as the re- 
maining spirals may simply have no current cool gas inflow. 
In this case, one would expect UV and Hi morphology to 
follow each other reasonably closely for the XUV disks but 
not for many of the others, as the star- formation in the inner 
disk is more closely related to the molecular phase (Bigiel 
et al. 2008). However, the existence of xuv disks pose an 
intriguing alternate possibility for the origin of the atomic 
hydrogen disk. Instead of primordial gas accreting onto the 
disk, the Hi disk could also be the byproduct of photodis- 
sociation of molecular hydrogen on the 'skins' of molecular 
clouds by the ultraviolet flux of the young stars in the XUV 
disk (see Allen et al. 1997; Allen 2002; Allen et al. 2004). 
This explanation has been explored in the stellar disks of 
several nearby galaxies (Heiner et al. 2008a, b, 2009, 2010). 
Gil de Paz et al. (2007b) calculate the time-scales (molec- 
ular gas dissociation and re-formation) involved but these 
are inconclusive regarding the origin of the xuv disk. In this 
scenario, one would expect the UV and H I morphologies to 
follow each other closely in all cases; in the outer disk, the 
UV flux from a few young stars would reach out to large 
areas of the low-column density gas to dissociate enough 
hydrogen to form the outer H I disk. Thus, the low- flux and 



low H I column density morphology -those defining the lim- 
its and extent of the XUV and H I disks- should show a close 
relation in parameters such as concentration, Gini, M20 and 
the effective radius. A third explanation is in terms of re- 
cent tidal interaction. A major merger often pulls gas out of 
the planes of galaxies and triggers star- forming events. Some 
anecdotal evidence (e.g., UGC 04862 is a late-stage major 
merger Torres-Flores et al. 2012) does point to this possible 
origin. In this scenario, one would expect that most H I disks 
hosting an XUV disk would be seriously tidally disrupted. 

In this series of papers, we have explored the quanti- 
fied morphology of available H i maps with the common pa- 
rameters for visible morphology; concentration-asymmetry- 
smoothness, Gini and M20 and Gm- In Holwerda et al. 
(2011c), we compare the Hi morphology to other wave- 
lengths, noting that the Hi and ultraviolet morphologies 
are closely related. In subsequent papers of the series, we 
use the H l morphology to identify mergers (Holwerda et al. 
201 Id), their visibility time (Holwerda et al. 2011a) and sub- 
sequently infer a merger rate from WHISP (Holwerda et al. 
2011b) as well as identify phenomena unique to cluster mem- 
bers (Holwerda et al. 201 le). 

In this paper, we explore the morphological link be- 
tween the Hi and xuv disks in the Westerbork Hi Survey 
Project (WHISP, van der Hulst et al. 2001; van der Hulst 
2002), a survey of several hundred H i observations of nearby 
galaxies. We complement this data with GALEX images to 
explore the morphology and typical scales of these maps. A 
direct and quantified comparison between the gas and ongo- 
ing star-formation morphology could help answer open ques- 
tions regarding the origin and nature of xuv disks: how do 
their respective sizes relate? Are their morphologies closely 
related in every case? Do their respective morphologies point 
to a dominant formation mechanism; gas accretion, pho- 
todissociation or tidal? Are xuv disks in morphologically 
distinct or typical H i disks? Are xuv disks embedded in H i 
disks that appear to be in an active interaction? What is 
the relation between UV flux and H I column density in the 
XUV disks, especially the outer disk? 

The paper is organized as follows; §2 gives the defini- 
tions of the quantified morphology parameters we employ, 
and §3 describes the origin the data. We describe the appli- 
cation of the morphological parameters and the results in §4 
and 5, and we discuss them in §6. We list our conclusions 
and discuss possible future work in §7. 



2 MORPHOLOGICAL PARAMETERS 

In this series of papers, we use the Concentration- 
Asymmetry-Smoothness parameters (CAS, Conselice 2003), 
combined with the Gini-M2o parameters from Lotz et al. 
(2004) and one addition of our own Gm- We have discussed 
the definitions of these parameters in the previous papers, 
as well as how we estimate uncertainties for each. Here, we 
will give a brief overview but for details we refer the reader 
to Holwerda et al. (2011c,d). 

We select pixels in an image as belonging to the galaxy 
based on the outer H I contour and adopt the position from 
the 2MASS catalog (Kleinmann et al. 1994) as the central 
position of the galaxy. Given a set of n pixels in each object, 
iterating over pixel i with value U , position Xi , yi with the 
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centre of the object at Xc,yc these parameters are defined 
as: 



C = 5 log (rso 7^20 ), 



(1) 



with r/ as the radial aperture, centered on Xc,yc containing 
percentage / of the light of the galaxy (see definitions of r/ 
in Bertin & Arnouts 1996; Holwerda 2005). ^ We include 
the rso, or "effective radius" in our catalog as well. 



A = 



— /i8o| 



(2) 



where /iso is the pixel at position i in the galaxy's image, 
after it was rotated 180° around the centre of the galaxy. 



S = 



(3) 



where Is is pixel i in a smoothed image. The type of smooth- 
ing (e.g., boxcar or Gaussian) has changed over the years. 
We chose a fixed 5" Gaussian smoothing kernel for simplic- 
ity. 

The Gini coefficient is defined as: 

G= ^ M2i-n-l)Iu (4) 
In{n — 1) 

where the list of n pixels was first ordered according to value 
and / is the mean pixel value in the image. 



M20 



(1^) ' for ^ili < 0.2/tot, 
V Mtot / 



(5) 



where Mi is the second order moment of pixel z; Mi = U x 
[{x — XcY + (y — UcY]- Mtot is the second order moment 
summed over all pixels in the object and M20 is the relative 
contribution of the brightest 20% of the pixels in the object. 
Instead of using the intensity of pixel z, the Gini parameter 
can be defined using the second order moment: 



Gm = 



Mn(n - 1 



-T.i{2i-n - l)Mi, 



(6) 



These parameters trace different structural characteris- 
tics of a galaxy's image but these do not span an orthogonal 
parameter space (see also the discussion in Scarlata et al. 
2007). 



3 DATA 

We use three data sets for this paper; H I radio observations 
from WHISP or THINGS, and ultraviolet observations from 
GALEX space telescope observations. 



3.1 WHISP Sample 

The starting dataset here is the 266 observations done as 
part of the Westerbork observations of neutral Hydrogen in 
Irregular and SPiral galaxies ( WHISP, van der Hulst et al. 
2001; van der Hulst 2002, ; 339 individual galaxies) that also 
have reasonable quality FUV and NUV data. The WHISP 



^ We must note that the earlier version of our code contained 
an error, artificially inflating the concentration values. A check 
revealed this to be Cnew = 0.38Coidi cind we adopt the new, 
correct values in this paper. 



observation targets were selected from the Uppsala General 
Catalogue of Galaxies (Nilson 1973), with blue major axis 
diameters > 2.^0, declination (B1950) > 20 degrees and ffux 
densities at 21-cm larger than 100 mJy, later lowered to 20 
mJy. Observation times were typically 12 hours of integra- 
tion. The galaxies satisfying these selection criteria generally 
have redshifts less than 20000 km/s {z < 0.07). WHISP has 
mapped the distribution and velocity structure of H I in sev- 
eral hundreds of nearby galaxies, increasing the number of 
Hi observations of galaxies by an order of magnitude. The 
WHISP project provides a uniform database of data-cubes, 
zeroth-order and velocity maps. Its focus has been on the 
structure of the dark matter halo as a function of Hubble 
type, the Tully-Fisher relation and the dark matter content 
of dwarf galaxies (Swaters et al. 2002b; Swaters & Balcells 
2002; Noordermeer et al. 2005b; Zwaan et al. 2005). Un- 
til the large all-sky surveys with new instruments are com- 
pleted, WHISP is the largest, publicly available data-set of 
resolved H I observations^ . We compiled a catalogue of basic 
data (position, radial velocity etc.) (see for details Holwerda 
et al. 2011b). 

3.1.1 Hi Data 

We use the highest available resolution zero-moment maps 
(beam size of -12" x 12"/sin(^)), from the WOW website 
and converted these to Mq/pc^ column density maps with 
J2000 coordinates and of the same size as the GALEX postage 
stamp (see for details Holwerda et al. 2011b). 



3.1.2 GALEX data 

To complement the Hi column density maps, we retrieved 
GALEX (Martin et al. 2005) postage stamps from http: 
//skyview. gsfc.nasa.gov, near- and far-ultraviolet (1350- 
1750 and 1750-2750 A). The majority of observations were 
taken as part of the AIS, MIS, NGS, etc. imaging surveys 
and hence vary in depth and signal-to- noise. Spatial resolu- 
tion is 4-6" , depending on position in the field-of-view. 

Of the 339 galaxies in WHISP, 266 had rehable FUV 
information and this cross-section is what we use for the 
remainder of the paper. The FUV images with H I contours 
overlaid are shown in Figures 1, 2, and Figures ?? ?? and 
all galaxies in appendix C {electronic edition). 



3.1.3 xuv classifications 

In the WHISP sample, Thilker et al. (2007) classified 22 
galaxies for the presence and type of XUV disks. The Type 1 
xuv disks are characterized by discrete regions out to high 
radii and Type 2 by an anomalously large, UV-bright outer 
zone of low optical surface brightness. Table 1 lists the XUV 
classifications for the WHISP sample; 22 from Thilker et al. 
(2007) and a single one from Lemonias et al. (2011). The 
latter study's lack of overlap is unsurprising as they focus on 
a more distant sample in complement to the Thilker et al.'s 
nearby sample {D < 40Mpc). We designate an XUV disk with 
both a Type 1 and Type 2 as Type 1/2 to discriminate in 



2 Available at "Westerbork On the Web" (WOW) project at AS- 
TRON (http : //www . astron . nl/wow/) . 
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the plots below and mark those galaxies in WHISP surveyed 
by Thilker et al. that do not have an XUV disk as Type 0. 

Figures 1 and 2 show the overlap between WHISP and 
the XUV disks identified by Thilker. The Type 1 collection 
includes a single merger remnant (UGC 4862) and some H I 
disks that appear slightly lopsided, tentatively supporting 
a tidal origin as one path to generate Type 1 disks. Figure 
3 shows the H i contours around those WHISP galaxies in- 
spected by Thilker et al. but not found to contain an XUV 
disk. There are two instances of clear ongoing mergers as 
well as a variety of Hi morphologies. 



3.2 THINGS 

Since there overlap between WHISP sample and the Thilker 
et al XUV classification sample is small, we also compare the 
THINGS survey's H l and FUV morphologies to their clas- 
sifications of XUV disks (Table 2). We use the morphological 
information from Holwerda et al. (2011c), with the concen- 
tration definition appropriately amended. The morphologi- 
cal classifications are based on the public datasets from the 
Hi Nearby Galaxy Survey (THINGS, Walter et al. 2008)^ 
and GALEX Nearby Galaxy Atlas (NGA, Gil de Paz et al. 
2007a)^, retrieved from MAST (http://galex.stsci.edu). 
We use the robust-weighted (RO) THINGS column density 
maps as these are similar in resolution as the native resolu- 
tion of GALEX for the contour of 0.3 x 10^° cm~^, approx- 
imately the spatial extent of the Hi disk. Table 2 lists the 
THINGS galaxies that have XUV disk classifications from 
Thilker et al. (2007). The benefits of the THINGS sample 
are that there is an XUV classification for all galaxies, the 
Hi and UV data are of comparable spatial resolution, and 
the GALEX observations are of uniform depth. Drawbacks 
are that is is an equally small sample as the WHISPoyerldip 
and there is only a single ongoing merger in this sample as 
it was constructed to be quiescent. 



4 APPLICATION OF QUANTIFIED 

MORPHOLOGY TO H I AND GALEX DATA 

Before we compute the quantified morphology parameters 
described above, we preprocess the images as follows. First, 
we cut out a postage stamp of both the H l maps and GALEX 
data centered on the UGC position. Similar to Holwerda 
et al. (2011b), we use a stamp size of 22 x 22 arcmin around 
this position. We then smooth the GALEX data to approxi- 
mately the same spatial resolution as the WHISP HI data (~ 
12"). We use the Hi map to define the part of the data over 
which the morphological parameters are computed over in 
the same manner as in Holwerda et al. (201 Id, b); a threshold 
of 10^° atoms/cm^ 

We then compute the above quantified morphology pa- 
rameters for the Hi map and far-, and near-ultraviolet im- 
ages (summarized in Tables ??, ?? and ??, full versions in 
the electronic edition). In addition to the above parameters, 
we compute the effective or half-light radius (-R50) for each 
image as an indication of disk size on the sky. 



Table 1. xuv disks classified by Thilker et al. (2007) and one from 
Lemonias et al. (2011) for the WHISP galaxies in their respective 
samples. 



Name 


UGC 


Type 


- 


1913 





- 


4165 





- 


4274 





- 


4273 





NGC 2541 


4284 


2 


- 


4325 





- 


4499 





NGC 2782 


4862 


1 


- 


5556 





NGC 3319 


5789 


2 


NGC 3344 


5840 


1 




6856 







7323 





NGC 4258 


7353 


1 




7651 





NGC 4559 


7766 


1 




7853 





NGC 4625 


7861 


1 


NGC 5832 


9649 


1/2 




10445 


1/2 




10791 


1 




12754 






Table 2. xuv disks classified by Thilker et al. (2007) for those 
THINGS galaxies in their sample. 



Name 


Type 


Reference 


NGC 628 


1 


Lelievre &; Roy (2000) 
Alberts et al. (2011) 


NGC 925 







NGC 2403 


1/2 




Holmberg-II 







M81A 







DD053 







NGC 2841 


1 


Alberts et al. (2011) 


NGC 2903 







HolmbergI 







NGC 2976 







NGC 3031 


1 




NGC 3184 







NGC 3198 


1 




IC2574 


2 




NGC 3351 







NGC 3521 







NGC 3621 


1 


Alberts et al. (2011) 


NGC 3627 







NGC 4736 





Zaritsky & Christlein (2007) 


DD0154 







NGC 4826 







NGC 5055 


1 


Alberts et al. (2011) 


NGC 5194 







NGC 5457 


1 




NGC 6946 







NGC 7331 







NGC 7793 








^ http : //www . mpia-hd . mpg . de/THINGS/ 
^ http://galex.stsci.edu/GR4/ 
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Figure 1. The few examples of xuv disks, Type 1 (individual UV complexes) as classified by Thilker et al. (2007) with complementary 
Hi data in WHISP. The tidal feature in UGC 04862 is studied in detail by Torres-Flores et al. (2012), who report high metallicities and 
1-11 Myr ages for the UV complexes. They speculate that their origin is either from the galaxies that formed UGC 04862 or that the 
current star-formation is a later generation, enriched system. 



5 ANALYSIS 

To explore in what type of H i disk xuv disks predominantly 
occur, we first compare the relative effective radii of the H i 
and UV disks, secondly compare the individual morphologi- 
cal parameters in both H I and UV against one another, and 
thirdly explore the UV and Hi morphologies alone, to ex- 
plore if xuv can be identified from UV morphology and to 
see if XUV disks are predominantly in merging or interacting 
disks. 



5.1 Effective Radii 

First we compare the relative sizes of the NUV and FUV 
disks and the H I disk t brought their effective radii, the ra- 
dius which contains 50% of the total flux. A naive expecta- 
tion would be that in the case of extended UV disks, these 
would be similar to the H I values. Figure 4 shows the ratio 
of the effective radii {Reff) for Hi and FUV or NUV. Nei- 



ther ratio singles out those identified by Thilker et al. (2007) 
as extended UV disks. Thus the ratio of effective radii is not 
a good metric to automatically select XUV disks. 

5.2 UVvs Hi Morphologies 

Figures 5 and 6 show the Concentration, Asymmetry, 
Smoothness, Gini, M20 and the Gm of the NUV and FUV 
compared to their values in the Hi maps, as computed over 
the same area and convolved to the same spatial resolution. 
We note that our approach is different from that of Lemonias 
et al. (2011) who use the optica/ concentration to identify the 
approximate galaxy type. 

Concentration of the UV and the Hi shows two dis- 
tinct sets of outliers. One where the UV concentration is 
more than three; these are relatively extended disks. Some 
of the XUV disks from Thilker et al. (2007) can be found 
here, but by no means all; two type 1/2 and one type 1 is 
not. However, six of the non-XUV disks are also above the 
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Figure 2. The two known examples of Type 2 xuv disks, a ring of UV outside the optical disk (top) and the two examples of Type 
1/2, a combination of Types 1 and 2, as classified by Thilker et al. (2007), which are also in the WHISP sample. The Hi morphology is 
remarkably flocculant and lacks grand spiral structure. 



Cfuv > 4 line. Figures 7 and 8 show that most Type 1 
in the THINGS sample are here. The second set of outliers 
have H l concentration indices exceeding Chi > 5. However, 
these are not necessarily the classical H I disks which extend 
well outside the stellar disk. High concentration values such 
as these often indicate recent tidal activity (Holwerda et al. 
201 Id). We will use the Concentration indices in FUV and 
H I to set our own supplementary definition of an extended 
disk; XUVc disks are defined by Cfuv > 4 and extended H I 
disks (xHi) by Chi > 5 in the other plots. We do so to iden- 
tify these outliers in the other plots to explore their nature, 
not to supplant the original detailed classification by Thilker 
et. al. The XUVc disks have UV flux distributed throughout 
their otherwise unremarkable H I disk and XHI disks are ex- 
tended in H I but not special in UV. If we compare the ratio 
of concentration indices in Figure 9 for the WHISP galax- 
ies, both the XUVc and XHI disks -unsurprisingly- stand out 
very clearly. This is confirmed by the ratios of THINGS con- 
centration index in Figure 10. We note, however, that if one 
applies these criteria blindly to the THINGS data (dashed 



lines in Figure 8), one would conclude that most THINGS 
galaxies are xuVc disks. Morphological criteria such as these 
need to be recalibrated when applied to different quality 
data. 

xuv and xuVc disks have much lower values of the M20 
parameter in the UV than the general WHISP population. 
There is a general trend between M20 and concentration 
(specifically for the 3.6 fiui Spitzer images Munoz-Mateos 
et al. 2009, Holwerda et al. in preparation), so extreme val- 
ues for this parameter are not surprising as such, especially 
if identified with a concentration criterion (see also Figure 
11. The lower value of M20 implies that regions that make 
up the XUV disk do not contribute to the brightest 20% of 
all the pixels in these UV images.; the general UV flux is 
more extended and hence the second order moment is less 
concentrated in the brightest few regions. The H I M20 val- 
ues however are quite within the normal range of WHISP 
galaxies. Thus, perhaps the UV M20 values offer an alterna- 
tive automated definition of XUV disks in combination with 
concentration. 
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The XUV disks show a median value for Asymmetry in 
FUV, but the full range in Hi. 

In UV smoothness, the XUV disks do not stand well 
apart from the general WHISP population. In the case of the 
G and Gm parameters, XUV disks avoid extreme values for 
the UVand generally fall below the Gm(IIi) =0.6 criterion 
for interactions. In the case of most of the H I morphological 
parameters, the values for the XUV disks are typical for the 
bulk of WHISP galaxies. 



Based on their relative Hi and UV morphologies, the 
previously identified or automatically identified XUV disks 
mostly stand out in their UV concentration and M20 values. 
In the remaining UV parameters and all Hi morphological 
values, they appear to be mixed in with the bulk of the 
WHISP (Figures 5 and 6) and THINGS (Figures 7 and 8). 
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Figure 4. The ratio of effective radii (i?5o) of the Hi over UV as a function of Hubble type as a function of type. 



Table 3. xuv disk selection criteria and their success (percentage of all bona-fide xuv disks included) and contamination rate (percentage 
of included objects that are non-xuv): 



Criterion 



WHISP 



THINGS 



XUV (%) Cont (%) Nr XUV (%) Cont (%) 

Cfuv > 4 30. 57 31 25 71 

|A- (-0.26667 *M2o) +0.12)/2| < 0.14 80 55 109 - 

|A - (-0.26667 *M2o)+ 0.12)1 < 0.14 - - - 87 56 
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Figure 9. The ratio of Concentration in Hi over that in UV as 
a function of Hubble type. In this case (as expected), the XUV 
and XHI disks stand out. 



Figure 10. The ratio of concentration of the Hi over UV (NUV 
and FUV) for the THINGS galaxies 



5.3 xuv disk UV morphology 

In this section we look at the UV morphological parameters 
in detail to explore if bona-fide xuv disks can be identified 
from their UV morphology alone -as computed over an area 
defined by the H i disk. For this we can use the 22 galaxies in 
WHISP and the THINGS sample, both classified by Thilker 
et al. 



Figure 11 shows the FUV parameters for the xuv disks 
identified by Thilker et al. (2007) as well as those disks in- 
spected but without an extended UV component for both 
samples. 

The above Concentration-only criterion (xuVc, Cfuv > 
4) pre-selects a number of the xuv disks but it would miss 
several bona-fide xuv disks as well as include many false 
positives. Using the information in two (or possibly more) 
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Figure 5. The morphological parameters of the Hi disk and NUV data; concentration (C), asymmetry (A), smoothness (S), the Gini 
index (G), and the second order moment of the brightest 20% of the pixels (M20), and the gini index of the second order moment (Gm)- 



parameters is how we identify mergers and may identify 
(candidate) XUV disks On inspection of the distributions of 
the NUV and FUV parameters (Figures ?? and ??), there 
are several combinations of parameters that offer the pos- 
sibility of pre-selecting xuv disks without too much loss: 
Gini- Asymmetry, M20- Asymmetry and Gini-CM. The goal 



is not to just select only the bona- fide, previously identi- 
fied xuv disks but also to exclude as many disks, that were 
previously classified as non-xuv as practical. This latter cri- 
terion immediately excludes Gini- Asymmetry and Gini-GM 
combinations. The combination of M20 and Asymmetry in 
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Figure 6. The morphological parameters of the Hi disk and FUV data; concentration (C), asymmetry (A), smoothness (S), the Gini 
index (G), and the second order moment of the brightest 20% of the pixels (M20), and the gini index of the second order moment (Gm)- 



FUV is promising and we define the following criterion to 
identify xuv disks in an FUV survey (6" resolution): 

\A - (-0.27 X M20 + 0.12)1 < 0.14, (7) 

which works well for the THINGS sample and we divide the 



linear relation by a factor two for the WHISP survey (to 
account for the ^-^12" resolution): 

\A - (-0.27 X M20 + 0.12)/2| < 0.14. (8) 

Both are illustrated in Figure 11 with dashed lines. Table 
3 lists the two criteria the simple concentration criterion 
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Figure 7. The morphological parameters of the H I disk and NUV for the THINGS galaxies with the xuv disk classification from Thilker 
et al. (2007). 



and the M20- Asymmetry one (adjusted for each survey) and 
their success rate (percentage of all bona-fide xuv disks in- 
cluded) and contamination rate (percentage of all classified 
objects that are non-xuv). The M2o-Asymmetry criterion 
selects most xuv but with a contamination of some 55% of 
non-xuv disks. A future search for xuv disks could therefore 
use this as a first cut before visual classification. Given the 
109 galaxies selected by this criterion (Table ??), 61 

of these should be bona-fide xuv disks, 23% of our sample 
of 266. 



5.4 Are xuv disks in interacting Hi disks? 

Figure 12 shows five of the morphological criteria we used 
or defined in Holwerda et al. (2011a) to identify those disks 
that are currently undergoing an interaction (see Figure ?? 
for a similar plot to Figure 5 in Holwerda et al. (201 Id) for 
a direct comparison). We mark those galaxies identified by 
Thilker et al. (2007) and those identified as extended UV or 
H I disks by the concentration parameters. 

Based on the H i morphologies, the Gm and C — M20 
criteria, those Hi disks that host an xuv disk are not in- 
teracting in the majority of cases. In the case of the other 



morphological criteria {A — G, A — M20, G — M20), this dis- 
tinction is not as clear-cut but xuv disks -either selected by 
Thilker et al or by their concentration- appear to reside in 
morphologically very typical H I disks. 

At the same time, extended H i disks are often interact- 
ing. This is in part because of our definition of an interac- 
tion is based on concentration but even in the Gm criterion, 
these are conspicuously separated from the bulk of the Hi 
disks. Typically, the bona-fide xuv disks selected by an Hi 
morphological selection of an ongoing merger appear to be 
Type 1, consistent with the first impression from Figures 3 
- 2. Thus, mergers appear to be one avenue to generate a 
Type 1 xuv disk. 

If we plot these -as a sanity check- for the THINGS 
sample, the xuv disks identified by Thilker et al. (2007) 
equally appear not to be interacting based on the 
Cm criterion, but several are based on G-A or M20-G crite- 
ria. Several of the non-xuv disks are interactions in A-M20. 
These criteria would need to be adjusted for the higher reso- 
lution of the THINGS sample (difficult to do without many 
mergers in the THINGS sample, as we noted in Holwerda 
et al. 2011c). 
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Figure 8. The morphological parameters of the Hi disk and FUV for the THINGS galdocies with the xuv disk classification from Thilker 
et al. (2007). Dashed lines are the Hi and UV concentration criteria from Figure 6. 



6 DISCUSSION 

Our first motivation for this study was to explore the possi- 
bility if xuv disks could be reliably found using the morpho- 
logical parameters in UV and H i, especially when computed 
over a wide aperture such as the Hi disk. Based on their 
morphology parameters in UV and Hi, XUV disks identified 
in Thilker et al. (2007) rarely stand out from the bulk of 
the WHISP sample, with the exception of M20-A and con- 
centration in UV. In the FUV, the M20-A relation appears 
reasonably successful in identifying xuv disks and the rate 
in the WHISP sample (23%) is consistent with the Thilker 
et al. (2007) and Lemonias et al. (2011) result, for a survey of 
late-types such as WHISP. The FUV values of Asymmetry 
and M20 could therefore be used to identify candidate xuv 
disks in the local Universe with GALEX or similar quality UV 
data or in deep rest-frame UV observations with HST. How- 
ever, these will have to be computed over a much larger area 
than commonly used (e.g., the de Vaucouleur diameter D25 
or the Petrosian radius) , such as the one defined by the outer 
H I contour. The quantified H i and UV morphologies do not 
correlate as closely as they did for the THINGS sample in 
Holwerda et al. (2011c) (their Figure 6). In our opinion, this 



is due to the wider range of environments probed by the 
WHISP survey as well as its lower resolution compared to 
THINGS. A closer relation may return for the inner stellar 
disk and possibly at higher spatial resolution. 

Our second motivation for this work on XUV disks was 
to explore their origin. Are they the (by)product of (a) cold 
flow accretion, (b) a gravitational interaction or mergers 
with another galaxy or (c) is photo-dissociation of molec- 
ular clouds by the UV regions of the xuv disks responsible 
for the Hi disks? 

The fact that the majority of the xuv disks do not reside 
in morphologically distinct Hi disks (Figure 6 and 8), pre- 
cludes the notion that most XUV disks originate from major 
interaction between galaxies. Few are selected by the pre- 
viously deflned interaction criteria for Hi disks. There are, 
however, some obvious exceptions (Figure 1). So our naive 
interpretation for the origin of Type 1 disks in Figure 1 is 
only partly supported by H I morphology and only holds for 
Type 1. 

Since extended UV (xuVc) and extended Hi disk (xHi) 
populations -as defined by their concentration parameter 
computed over the same area- do not coincide (Figure 6, 
top left panel), a direct link, such as the photo-dissociation 
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Figure 11. Concentration, Asymmetry and M20 in FUV for the WHISP (left) and THINGS (right) samples with the Thilker et al. 
classifications marked. Cpuv > 4 does not select cleanly the xuv disks previously identified by Thilker et al. The M2o-Asymmetry 
criteria (green dashed lines in the middle panels), do select a relatively clean sample of xuv disks from both WHISP and THINGS. 



origin of the Hi disk, appears less likely. A closer UV-Hi 
relation in morphology may well still hold within the inner 
stellar disk, especially since the photo-dissociation scenario 
was developed to explain the higher column densities of H i 
in this environment (Allen et al. 1997, 2004; Heiner et al. 



2008a,b, 2009, 2010). However, the existence of both the 
Cfuv > 4 disks without complementary extended Hi disks 
as well as extended (Chi > 4) Hi disks without a comple- 
mentary xuv disk appears to be contradictory to the photo- 
dissociation origin of XHI disks. Two considerations are still 
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Figure 12. The morphological selection criteria (dashed lines) for Hi applied to the WHISP sample. Disks are interactions above and 
to the right of the dashed lines, xuv disks identified either by Thilker et al. or high FUV or H I concentration are marked. 



in its favor however, Gogarten et al. (2009) showed that UV 
emission does not fully trace the lower-lever star-formation 
in the low-column density environments of the outer disk. 
Thus, the ionizing stars may be lower-mass Main Sequence 
stars and simply remain undetected in GALEX imaging due 
to their extreme low surface-brightness. And with higher 



resolution and more sensitive H i and UV imaging, a closer 
relation may still return in their quantified morphology (wit- 
ness the THINGS sample). Thus, we cannot completely rule 
out a photo-dissociation scenario explaining the relation be- 
tween Hi and UV disks. Similarly, Gil de Paz et al. (2007b) 
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find their time-scale argument for a photo-dissociation ori- 
gin of XUV disks to be inconclusive. 

This leaves cold flows as a likely origin for most xuv 
disks. A number of these Type 1 xuv disks may be the 
result of a later-stage full merger such as UGC 04862 but 
judging from the Hi morphologies, this is small fraction of 



the full xuv or xuVc sample. Cold flows occur as both a 
steady trickle along a fllament of the Cosmic Web as well 
as discrete accretion in the form of the cold gas originally 
associated with a small satellite. 

There are indirect signs that cold flows are indeed the 
most common cause for XUV disks. For example, at least 
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three of the nine THINGS galaxies with an XUV disk (M83, 
NGC 4736 and NGC 5055) are currently cannibalizing a 
small companion (Martmez-Delgado et al. 2010). Satellites 
embedded in the gas filament are one of two fueling mech- 
anisms of cold flow; discrete and trickle cold flows (Keres 
et al. 2005). The acquired cold gas accompanying the can- 
nibalized companion is a likely trigger for the formation UV 
complexes. Similarly, those disks selected by their extended 
UV emission (Figure ??), often appear to have flocculent, 
low-column density Hi morphologies, which may point to 
recent cold accretion. A few cases (UGC 7081, 7651, and 
7853) are obviously interacting major mergers with accom- 
panying H l/UV tidal features. In contrast, the XHI disks are 
nearly always small galaxies with an obvious close compan- 
ion (Figure ??), which makes a tidal origin likely for their 
very extended H i disk. 

The quantified morphologies of the WHISP H i and UV 
disks are not conclusive evidence that cold flows are respon- 
sible for the XUV complexes but certainly suggest this is the 
most likely origin for most, with a few exception due to ma- 
jor interactions. 



7 CONCLUSIONS 

Based on the quantified morphology of the H i and UV maps 
of the WHISP sample of galaxies, we conclude the following: 

1. There are distinct galaxy populations that stand out 
by their high concentration values in FUV or H I (xuVc and 
XHI thoughout the paper). These population do not overlap. 
Some known XUV disks are in the XUVc sample (Figures 6-a 
6-e, 8-a and 8-c).To compute this concentration index, the 
outer Hi contour is needed to delineate the extent of the 
disk. 

2. The fact that relative dilute UV disks {Cfuv > 4) and 
H I disks {Chi > 5) population do not overlap at all and the 
lack of close morphological relations in any other parameters 
together suggest that a common small-scale origin for the 
UV and H I disks, such as a phot odissociat ion of molecular 
hydrogen scenario, is less likely for XUV disks (but may still 
very much hold for the inner disks). 

3. Asymmetry and M20 can be used in combination to se- 
lect XUV disks with reasonable reliability (80% included) but 
substantial contamination (55%), when properly calibrated 
for the survey's spatial resolution (equation 7 and 8) and 
the morphologies computed over a large enough aperture. 
With this selection, one can find the number of XUV disks in 
a survey or candidates for visual classification (Figure 11). 

4. Based on the morphology of the H I disk in which they 
occur, XUV disks appear not to occur often in tidally dis- 
turbed gas disks (Figures 12 and 13). 

5. In a few cases, the XUV disk is the product of a major 
merger; for example UGC 04862, identified by Thilker et al. 
(2007) and UGC 7081, 7651, and 7853, identified by their 
large FUV concentration. This appears to be an avenue to 
form some of the Type 1 XUV disks. 

6. The Hi morphology and anecdotal evidence for small 
satellite cannibalism all point to a third mechanism for the 
origin of most XUV disks; cold fiow accretion. 

With the emergence of new and refurbished radio ob- 
servatories in preparation for the future Square Kilometre 



Array (SKA; Carilli & Rawlings 2004), a new window on the 
21 cm emission line of atomic hydrogen gas (Hi) is opening. 
The two SKA precursors, the South African Karoo Array 
Telescope (MeerKAT; Booth et al. 2009; Jonas 2007; de Blok 
et al. 2009), and the Austrahan SKA Pathfinder (ASKAP; 
Johnston 2007; Johnston et al. 2007, 2008a,b, 2009) stand 
poised to observe a large number of Southern Hemisphere 
galaxies in Hi in the nearby Universe (z<0.2). In addition, 
the Extended Very Large Array (EVLA; Napier 2006) and 
the APERture Tile In Focus instrument (APERTIF; Ver- 
heijen et al. 2008; Oosterloo et al. 2009) on the Westerbork 
Synthesis Radio Telescope (WSRT) will do the same for the 
Northern Hemisphere. The surveys conducted with these 
new and refurbished facilities will add new, high-resolution, 
H I observations on many thousands of galaxies. In the case 
of WALLABY (Koribalski et al. in preparation), these will 
be of similar quality spatial resolution as the WHISP survey. 

Combining these with existing UV observations from 
GALEX, we can explore the relation between the morphology 
of the atomic hydrogen and ultraviolet light for much greater 
samples and in much greater detail. An application of the 
Asymmetry- M20 identification of xuv disks in the GALEX 
Nearby Galaxy Atlas (Gil de Paz et al. 2007a) or a sample 
similar to Lemonias et al. (2011) could reveal additional ex- 
amples but the combination with deep H I observations can 
prove the link with cold flows for the majority of XUV disks. 
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